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ABSTRACT 

We present the intrinsic spectral energy distribution (SED) of the Narrow-Line Seyfert 1 galaxy (NLS1) 
Ark 564, constructed with contemporaneous data obtained during a multi-wavelength, multi-satellite observing 
campaign in 2000 and 2001. We compare this SED with that of the NLS1 Ton SI 80 and with those obtained for 
Broad-Line Seyfert Is to infer how the relative accretion rates vary among the Seyfert 1 population. Although the 
peak of the SED is not well constrained, in our parameterization most of the energy of this object is emitted in the 
10-100eV regime, constituting roughly half of the emitted energy in the optical/X-ray ranges. This is consistent 
with a primary spectral component peaking in the extreme UV/soft X-ray band, and with disk-corona models, 
hence high accretion rates. Indeed, we estimate that m w 1. We also address the issue of the energy budget in 
this source by examining the emission lines observed in its spectrum, and we constrain the physical properties 
of the line-emitting gas through photoionization modeling. The available data suggest that the line-emitting gas 
is characterized by logn w 11 and \ogU w 0, and is stratified around \ogU w 0. Our estimate of the radius of 
the H/3-emitting region R B ^ R ~ 10 ± 2 lt-days is consistent with the 7? B L R -luminosity relationships found for Sy 1 
galaxies, which indicates that the narrowness of the emission lines in this NLS1 is not due to the Broad-Line 
Region being relatively further away from the central mass than in BLSls of comparable luminosity. We also find 
evidence for super-solar metallicity in this NLS1. We show that the emission lines are not good diagnostics for 
the underlying SEDs and that the absorption line studies offer a far more powerful tool to determine the ionizing 
continuum of AGNs, especially if comparing the lower- and higher-ionization lines. 

Subject headings: galaxies: active - galaxies: individual (Arakelian 564) - galaxies: nuclei - galaxies: Seyfert - 
galaxies: NLS1 - galaxies: emission lines 



1. INTRODUCTION 

The population of Seyfert 1 galaxies has a widely used 
sub-classification into Narrow-Line Seyfert 1 galaxies (NLSls) 
and Broad-Line Seyfert 1 galaxies (BLSls). While this 
classification appears to make an arbitrary distinction based 
on the widths of the optical emission lines (NLSls having 
FWHM(H/3) < 2000 km s" 1 , Goodrich 1989), this is in fact an 
extremely useful scheme since the X-ray properties of the two 
subclasses are systematically different. As a class, NLSls show 
rapid and large-amplitude X-ray variability (Boiler, Brandt, & 
Fink 1996, hereafter BBF96; Turner et al. 1999b), with the ex- 
cess variance (Nandra et al. 1997a) typically an order of mag- 
nitude larger than that observed for samples of BLSls with 
the same luminosity distribution (Turner et al. 1999b; Leighly 
1999a). Analogously, the spectral properties also vary across 
the Seyfert population with NLSls showing systematically 
steeper spectra than those of BLSls in both the soft and hard 
X-ray bands (Boiler, Brandt, & Fink 1996; Brandt, Mathur, & 



Elvis 1997; Turner, George, & Nandra 1998; Leighly 1999b; 
Vaughanetal. 1999). 

One increasingly popular hypothesis to explain the differ- 
ences in X-ray properties across the Seyfert population is that 
NLSls have relatively low masses for the central black hole 
compared to BLSls with similar luminosities. Smaller black- 
hole masses naturally explain both the narrowness of the op- 
tical emission lines, which are generated in gas that has rela- 
tively small Keplerian velocities, and the extreme X-ray vari- 
ability, since the primary emission would originate in a smaller 
region around the central engine (e.g., Laor et al. 1997). Given 
that NLSls have comparable luminosity to that of the BLSls, 
Pounds, Done, & Osborne (1995) suggested that they must 
be emitting at higher fractions of their Eddington luminosity, 
hence higher fractional accretion rates (m = M /Mndd) are also 
required. The closer the luminosity is to the Eddington limit 
(and the lower the black-hole mass), the greater the fraction 
of the energy emitted by the accretion disk in the soft X-rays 
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(Ross, Fabian, & Mineshige 1992). Thus NLSls might be ex- 
pected to show disk components which peak at higher ener- 
gies than for BLSls. Pounds, Done, & Osborne (1995) and 
Maraschi & Haardt (1997) noted that soft photons from the disk 
may Compton-cool hard X-rays from the corona, and cause the 
observed steep spectra. 

Given the dependence of the peak energy of the disk emission 
on the accretion rate, the spectral energy distribution (SED) of 
an active galactic nucleus (AGN) will provide critical informa- 
tion about accretion rates and conditions close to the disk. In 
particular, examination of the SED of a NLS1, and comparison 
with that obtained for BLSls will offer insight into the relative 
accretion rates across the Seyfert population. However, measur- 
ing the SED requires observations taken simultaneously over a 
long wavelength baseline stretching from infrared to hard X- 
ray wavelengths, and the determination of the UV/X-ray con- 
tinuum in AGN is still difficult due to the severe attenuation by 
even small amounts of Galactic interstellar gas along the line- 
of-sight. Another complication comes from intrinsic reddening, 
i.e., reddening associated with the active nucleus itself. Indeed, 
the steeper UV/blue continua observed in NLSls (when com- 
pared to AGN spectrum composites) can be attributed at least 
in part to reddening, though the ionization of the absorbing ma- 
terial and its location with respect to the accretion source is still 
not well determined (Constantin & Shields 2003). 

Arakelian 564 (Ark 564, IRAS 22403+2927, MGC +05-53- 
012) is a bright, nearby NLS1 galaxy, with z = 0.02467 and 
V = 14.6 mag (de Vaucouleurs et al. 1991), and a mean 2- 
10 keV luminosity L2-10 kev w 2.4 x 10 43 ergs s" 1 with flux vari- 
ations of a factor of a few in a few thousand seconds (Turner et 
al. 2001, hereafter Paper I). It was the object of an intense mul- 
tiwavelength monitoring campaign that included simultaneous 
observations from ASCA (2000 June 1 to July 6, Paper I; Pounds 
et al. 2001; Edelson et al. 2002), XMM-Newton (2000 June 
17, Vignali et al. 2003), Chandra (2000 June 17, Matsumoto, 
Leighly, & Marshall 2002), HST (2000 May 9 to July 8, Collier 
et al. 2001, Paper II; Crenshaw et al. 2002, Paper IV), FUSE 
(2001 June 29-30, Romano et al. 2002, Paper V), and from 
many ground-based observatories as part of an AGN Watch 12 
project (1998 November to 2001 January, Shemmeret al. 2001, 
Paper III). Ark 564 has shown a strong associated UV absorber 
(Crenshaw et al. 1999, Paper II; Paper IV; Paper V). There are 
indications that it also possesses a warm X-ray absorber, as seen 
by the narrow absorption lines of O VII and O VIII detected in 
a Chandra spectrum (Matsumoto, Leighly, & Marshall 2002), 
and that the UV and X-ray absorbers in Arakelian 564 are phys- 
ically related, possibly identical, and may be spatially extended 
along the line of sight (Paper V). 

In this paper we present a contemporaneous SED of Ark 564, 
based on the extensive monitoring of 2000. In §2 we describe 
the observations and data reduction, summarize the main results 
of the monitoring campaign, and describe the adopted method 
for reddening correction. In §3 we present the SED of Ark 564. 
In §4 we constrain the mean physical properties of the line- 
emitting gas through photoionization modeling. In §5 we dis- 
cuss some implications of our investigation. Finally, our results 
are summarized in §6. 



2. OBSERVATIONS 

2.1. Data Reduction 

Table 1 summarizes the log of the observations of Ark 564 
used in this study. Column 1 and 2 report the observatory and 
instrument which obtained the data; Column 3 shows the date 
in which the data were obtained; Column 4 lists the ranges of 
energy and wavelength we used for our work; Column 5 lists 
the total exposure times and the slit size, along with other rele- 
vant notes for a particular observation; finally, Column 6 reports 
the references to the papers where the data were first published. 
Most of the data have been published already, with the excep- 
tion of the XMM data which will be published in Vignali et al. 
(2003), therefore here we will only briefly summarize their re- 
duction and analysis, referring the reader to the original papers 
for further details, and to §2.2 for a summary of the main re- 
sults. 

The ASCA data (Paper I) were obtained starting from 2000 
Jun 1 (Julian Date 2451697.024) for a total exposure time of 
<~ 2.98 Ms, and were reduced using standard techniques as de- 
scribed in Nandra et al. (1997a) with the methods and screening 
criteria utilized by the Tartarus 13 database (Turner et al. 1999b). 
Data screening yielded an effective exposure time of ~ 1.11 
Ms for the SISs and — 1.29 Ms for the GISs. The degradation 
of the low energy response of the SIS detectors was corrected 
for with the method of Yaqoob et al. (2000) 14 , i.e. we param- 
eterized the efficiency loss with a time-dependent absorption 
(Vh(SISO) = 7.5 x 10 20 cm" 2 , V H (SIS1) = 1.05 x 10 21 cm" 2 , Pa- 
per I). 

FUSE observed Ark 564 for 63 ks starting from 2001 June 
29 (Paper V). The observations consisted of 24 exposures per- 
formed in photon address (time-tag) mode through the 30" x 
30" low-resolution (LWRS) aperture. The data were reduced 
with CalFUSE (version 2.0.5) 15 as a single continuous expo- 
sure. As a result of a high-voltage anomaly and data screening, 
the effective on-source times were 41 ks in Detector 1A, 39 ks 
in Detector IB, 58 ks in Detector 2 A, and 62 ks in Detector 2B. 
We discarded from further analysis the SiCl A and SiClB spec- 
tra because of a high- voltage anomaly, and the LiFIB spectrum 
since it showed wavelength-dependent differences in flux of up 
to 30-50 % compared to the LiFl A. We rebinned the full spec- 
trum in a linear wavelength scale using 0.6A bins (100 pixels, 
effective resolution of 20 km s" 1 ). Comparison with the flux 
level of the HST spectrum indicates that in June 2001, when 
the FUSE spectrum was obtained, the source was ~ 1.3 times 
brighter than in May-July 2000, when the HST spectrum was 
obtained. We therefore scaled the FUSE fluxes by 0.75 to match 
the HST data. As noted in Paper V, this is not inconsistent with 
a combination of effects such as flux intercalibration uncertain- 
ties and, most importantly, source flux variability. 

The HST data (Paper II and IV) were obtained with the Space 
Telescope Imaging Spectrograph (STIS) in 46 visits between 
2000 May 9-July 8 (the first five visits were separated by inter- 
vals of 5 days, the remaining by 1 day). The spectra were ob- 
tained through the 52" x 0"5 slit and the low-resolution G140L 
and G230L gratings which yield a spectral resolution of <~ 1 .2 A 
and3.2Ain the 1 150-1730 A and 1570-3 150 A ranges, respec- 
tively. The data were reduced using the IDL software 



All publicly available data and complete references to published AGN Watch papers can be found at http://www.astronomy.ohio-state.edu/~agnwatch. 
http://tartarus.gsfc.nasa.gov. 

see http://lheawww.gsfc.nasa.gov/~yaqoob/ccd/nhparam.html. 
See http://fuse.pha.jhu.edu/analysis/calfuse.html. 
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FIG. 1.— Quasi-simultaneous SED of Ark 564 (§2.1). The data have not 
been corrected for reddening and are in the observed frame. 

developed at NASA's Goddard Space Flight Center for the STIS 
Instrument Definition Team (Lindler 1998). The spectra have 
been corrected for small wavelength intercalibration uncertain- 
ties following Korista et al. (1995). The uncertainty in the rel- 
ative wavelength calibration is on the order of 0.6 A and 1 .7 A 
for the G140L and G230L gratings, respectively. A separate 
mean spectrum was created for the G140L and G230L grating 
separately, given the different resolutions. 

In the optical we combined two spectra. The first one, which 
covers the 3 1 70-4 1 60 A wavelength range with a mean spectral 
resolution of 0.62 A, was obtained in 1980 at Lick Observatory 
(D. E. Osterbrock 2002, private communication). The second 
one is the mean of the spectra taken between 1998 Nov to 2001 
Jan at the Tel Aviv University Wise Observatory (Paper III, res- 
olution of ^10 A). The host galaxy starlight contribution has 
been estimated by measuring its flux through PSF fitting to field 
stars in V-band images of the galaxy taken at Wise Observa- 
tory, which corresponds to <~ 40 % of the total light at 5200 A, 
i.e. F g . di = 2.4 x 10" 15 ergs s" 1 cm" 2 A" 1 (Paper III). Given the 
limited resolution (dominated by a seeing disk of <~ 2. "5), it 
was not possible to separate the components of the host galaxy 
(bulge and bar) from PSF, hence we subtracted from the mean 
spectrum a constant host contribution of F gii \- No scaling was 
necessary between the two spectra, in agreement with the low- 
amplitude variations of the continuum found in Paper II (~ 6% 
over a month-long observation) and III (~ 10% during a two- 
year monitoring). Indeed, there is growing evidence that while 
NLS 1 s show very rapid and giant X-ray variability, they show 
only slow and minor optical variability (O. Shemmer et al 2003, 
in preparation). 

To extend the SED in the IR, we derived four continuum 
points between <~ 10000 and 24000 A from Rodriguez- Ardila 
et al. (2002a,b), which were obtained on 2000 Oct 11 and 13 
with the NASA 3m IRTF telescope and the SPEX spectrome- 



FlG. 2. — (a): FUV/Optical rest-frame spectrum of Ark 564, obtained 
combining the FUSE spectrum (June 2001; 1000-1175 A), the HST G140L 
and G230L mean spectra (May-July 2000; 1175-3 143 A), the Lick spec- 
tram (1980; 3170-4160 A), and the Wise mean spectrum (1998 Nov-2001 
Jan; 4160-7790 A). The FUSE spectrum has been scaled by a factor of 0.75 
to match the flux level of the HST spectrum; the Wise spectrum has been 
corrected for the galaxy contribution (F gii i = 2.4 X 10~ 15 ergs s _1 cirr 2 A -1 ; 
Shemmer et al. 2001). (b): Full spectrum after correction for reddening us- 
ing a standard Galactic extinction curve with E(B — V) = 0.03 mag plus the 
intrinsic extinction curve that Crenshaw et al. (2002) derive for Ark 564 and 
E(B-V) = 0.14 mag. The data have been corrected for redshift (z = 0.02467; 
de Vaucouleurs et al. 1991). The extinction-corrected spectrum (b) has been 
offset by 1 X 10~ 14 ergs s _1 cnT 2 A -1 for clarity. 

ter. We also retrieved archival IRAS flux measurements at 
12, 25, 60, and 100 fim (Moshir et al. 1990) through the 
NASA/IPAC Extragalactic Database (NED). 

Figure 1 shows the contemporaneous SED of Ark 564 before 
correction for intervening (and intrinsic, see §2.3) absorption is 
applied. We note that while the HST and Wise spectra are si- 
multaneous (as well as simultaneous with the ASCA spectrum), 
the FUSE and Lick spectra were obtained one year later, and 
20 years earlier, respectively. The FU V/optical rest-frame spec- 
trum of Ark 564 covering the 1000-7790 A wavelength region 
is presented in Figure 2 (labeled as (a)). 

2.2. Summary of results from the multi-waveband observations 

The continuum fit to the mean ASCA spectrum (with a power- 
law model modified by Galactic absorption, Nh = 6.4 x 10 20 
cm -2 , Dickey & Lockman 1990) yields a slope T = 2.538 ± 
0.005 (Paper I). The strong excess of emission observed be- 
low 2 keV was parameterized as a Gaussian of peak energy E = 
0.57 ± 0.02 keV and mean equivalent width (EW)= 1 10!j^ eV. 
The soft hump component is also found to be variable in flux 
on timescales as short as 1 day and in shape on timescales as 
short as a few days (Paper I). Parameterization of the soft ex- 
cess as a black-body yields a temperature 7 = 1.8 x 10 6 K and 
luminosity Lbb = 2.48 x 10 38 ergs s" 1 (Paper I). A strong, ion- 
ized (E ss 7keV) Fe Ka line is detected, which shows varia- 
tions in flux and EW on timescales as short as a week (Paper I). 



4 



Romano et al. 2002 



lOOAtm lOAtm ljum lOOOA 0.1 keV IkeV 



45 - 



44 



O 



43 



T 



J i L 



12 



13 



14 




16 
i/(Hz) 



19 



FIG. 3. — Reddening corrected, rest-frame SEDs of Ark 564 (see §3). The 
optical/FUV data are the spectrum labeled (b) in Figure 2. The X-ray data are 
drawn from Paper I, while the IR data are from IRAS and IRTF (§2.1). The 
short-dashed line shows the power-law fit of the optical/FUV data (spectral in- 
dex a = 0.42 ± 0.01) extrapolated to higher energies. The long-dash-dot line 
is the power-law fit to hard X-ray data (XPL, a ASCA = 1 .538 ± 0.005, §2.2) ex- 
trapolated to lower energies. The short-dash-dot line is the black-body model 
for the soft X-ray excess (BB, §2.2). The solid line is a more conservative 
spectral energy distribution (SEDB) that connects the FUV and the soft X-ray 
data with a simple power law {a = 1.08, §3). 

The FUSE spectrum is dominated by the strong emission 
in the O VI AA1032, 1038 resonance doublet, and heavily sat- 
urated absorption due to Hydrogen Lyman- Werner bands, O VI 
and C III are observed at velocities near the systemic redshift 
of Ark 564 (Paper V). The available data suggest that the UV 
and X-ray absorbers are physically related, possibly identical, 
and spatially extended along the line of sight, and character- 
ized by total column density and ionization parameter logiVn 
(cm -2 ) ss 21 and logi/ m -1.5. The absorbing gas is in a state 
of outflow with respect to the nucleus and carries out a kinetic 
luminosity about one order of magnitude smaller than the ob- 
served radiative luminosity of the source (Paper V). 

Variations in the UV continuum flux (1365 A, Paper II) are 
well correlated with variations in both the hard (2-10keV) and 
soft (0.75-2keV) X-ray fluxes and with the soft hump flux. 
No significant lags are detected between the variations in the 
X-ray and UV bands (Paper I). The variations of the contin- 
uum at 3000 A (Paper II) and at 4900 A (Paper III) lag behind 
those at 1365 A by ~ 1 day and 1.8 days, respectively. These 
UV/optical delays were interpreted as evidence for a stratified 
continuum reprocessing region, possibly an accretion disk. 

The variations of the Lya emission line, which lag the vari- 
ations of the continuum at 1365 A by 3 days, were used, in 
conjunction with the line width, to determine the virial mass 



of the central black hole, M < 8 x 10 6 M Q (Paper II). This esti- 
mate is uncertain due to the low amplitude of the Lya emission 
line variations (1%). However, the estimate in Paper II agrees 
with the one obtained by Pounds et al. (2001) based on a power 
spectrum analysis of X-ray variability. The black hole mass 
and 5100 A luminosity of Ark 564 are consistent with the hy- 
pothesis that NLSls have lower black hole masses and higher 
accretion rates than BLSls of comparable luminosity. The low 
level variability observed in the emission lines is also different 
from most Seyfert 1 galaxies, which characteristically display 
variations of 10 % on similar timescales. 



2.3. Reddening Correction in the Optical/FUV 

Given the indications (Paper IV) that strong intrinsic neu- 
tral absorption is present in Ark 564 in excess of the Galactic 
absorption, special care has been taken in correcting the data 
for reddening. We used a standard Galactic extinction curve 
with E(B-V) = 0.03 mag plus the intrinsic extinction curve 
that Crenshaw et al. (2002, Paper IV) derive for Ark 564 and 
E(B-V) = 0. 14 mag. The HST extinction correction was extrap- 
olated linearly into the FUSE band, as suggested by Hutchings 
& Giasson (2001) and Sasseen et al. (2002). The effect of red- 
dening correction in the optical/UV bands is presented in Fig- 
ure 2, where the observed spectrum in the 1000-7790 A wave- 
length region (labeled as (a)) is compared to the absorption- 
corrected one (labeled as (b)). 



16 The fit was made to 9 bands: A : 



3. THE INTRINSIC SED OF ARK 564 

Figure 3 shows the Ark 564 data in the IR/X-ray range. A 
power-law fit of the continuum in the optical/FUV region 16 
yields Fx oc A~ L58±001 , hence spectral index a = 0.42 ±0.01 
(specific flux /„ oc v~ a ). Extrapolation of this power law in 
the X-ray regime greatly overpredicts the X-ray flux (dashed 
line). Analogously, the hard X-ray continuum slope (ccasca = 
r-1 = 1.538 ±0.005, §2.2, long-dash-dot line) extrapolated to 
the lower energies overpredicts the optical/FUV flux, as previ- 
ously noted by Walter & Fink (1993). Clearly, both the opti- 
cal/FUV and the X-ray power-laws must break at some energy 
between the FUV and soft X-ray. With the adopted redden- 
ing correction (§2.3), the spectral energy distribution peaks at 
~50eV. 

Table 2 summarizes some relevant data from the spectral en- 
ergy distribution derived using the simple parameterization of 
the combination of an optical/FUV power law with a = 0.42 
breaking at <~ 50 eV to ocasca = 1-538 (hereon SED A), as de- 
scribed above, as well as the IRAS and IRTF data points (§2.1). 
Column (1) is the rest wavelength/energy, Column (2) and (3) 
list the observed and reddening corrected values of vL v , re- 
spectively. We estimate that the number of ionizing photons is 
g(SEDA) w 10 55 photons s" 1 . 

We also considered a more conservative spectral energy dis- 
tribution (the solid line in Figure 3, hereon SEDB) in which the 
the FUV and the soft X-ray data are connected with a simple 
power law (a = 1.08), i.e, the combination of an optical/FUV 
power law with a = 0.42 breaking at 1000A to a w 1, then again 
breaking at w 0.8 keV to (Xasca = 1 -538. The number of ionizing 
photons for SEDB is g(SEDB) w 5 x 10 54 photons s" 1 , which 

1005-1007, 1029.5-1030.5, 1101-1107, 1114-1118, 1155-1180, 1350-1380, 1460-1500, 1620-1660, and 7040-7050 A. The 



uncertainties are purely statistical. The continuum fit in Paper II was performed only on the HST data corrected for Galactic reddening (E(B- V) = 0.06 mag) and 
produced F x oc A -0 - 88 * 01 ; this is comparable to our value for this reddening, Fx oc A~° ■ 73 ± - 01 , s i nC e the uncertainties are underestimated in both cases. 
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is consistent with what Crenshaw et al. (2002) found 17 . 

4. PHOTOIONIZATION MODELING 

In this section, we address the issue of the energy budget 
of Ark 564 by examining the emission lines observed in its 
spectrum, and we deduce the mean physical properties of the 
line-emitting gas, Hydrogen density n and ionization param- 
eter U (U = Q/4-k r 2 n c, r being the distance to the ionizing 
source and c the speed of light) through photoionization mod- 
eling given the assumed spectral energy distributions. We ex- 
pect a range of ionization to exist throughout the BLR, to accu- 
rately describe which multi-zone modeling would be required; 
however, here we use a single density and ionization parameter 
modeling to derive the mean BLR properties. For each of the 
input continua (§4.2), we considered a total hydrogen density of 
n = 10 9 , 10 11 , 10 12 , and 10 13 cm" 3 , and calculated the predicted 
intensities of the major emission lines, for a range of ionization 
parameters between logf/ = -4 and logf/ = 2. In the case of 
"table agn"(§4.2), we specified a grid of n and U values. For 
the other input continua, we normalized the SEDs with respect 
to the measured X-ray luminosity in the absorption-corrected 
rest-frame 2-10keV energy range (Z-2-iokeV = 2.4 x 10 43 ergs 
s" 1 ; Paper I), and specified the radius of the cloud, thus ob- 
taining U . In the case of SEDA and SEDB, using an observed 
spectral energy distribution assumes that the gas responsible for 
the emission lines sees the same ionizing continuum as the ob- 
server does; therefore, we expect SEDA and SEDB to yield 
more realistic predictions of the emitted line spectrum. 

4. 1 . Emission Line Fluxes 

The fluxes (relative to H/3) of the most prominent emission 
lines in the 1 150-68 17A wavelength range have been published 
in Table 2 of Paper IV; here we report a selection of them in 
Table 3 (Column (2), relative to Lya). Because of the NLS1 
nature of Ark 564, the contribution from the BELR and NELR 
are strongly blended together and those line ratios include both 
components. The measured fluxes were corrected for reddening 
using the continuum reddening curve (Crenshaw et al. 2002), 
given the similar extinctions for the continuum and emission 
lines; the errors are propagated in quadrature from the ones 
listed in Paper IV, and they include photon noise, continuum 
placement errors, and reddening errors. 

Columns (2) and (3) of Table 3 also report the fluxes of 
C IIIA977 and the O VlAA1032,1038 doublet, which can help 
better constrain the value of n and U. These lines were mod- 
eled in Paper V, hence both a broad and a narrow component are 
available (denoted with BEL and NEL, respectively; Columns 
(2) and (3)). In order to compare with the other line ratios, we 
needed to account for the different contribution of BELR and 
NELR gas to H/3, which we estimated as follows. We assumed 
that the ratio of the narrow component of H/3 and [O III] A5007 
in NGC 5548, i.e., 0.12 ±0.01 (Kraemer et al. 1998), can be 
used for Ark 564; we scaled it to the observed [O III] A5007 flux 
for Ark 564, F([0 III]A5007) = (2.4±0.1) x 10 13 ergs s" 1 cm" 2 
(Paper III), obtaining the NEL ratio F(H/3)/F([0 III]A5007) in 
Ark 564. Thus, we estimated that roughly 75% of the total H/3 
flux is from the BEL, and bracketed this value between 50% and 
100%, given that the BEL and NEL contributions to C IIIA977 
and O VIAA1032, 1038 are the same, and that the total flux is the 

17 Note that SEDB corresponds to SED2 in Paper V. 

18 http://www.nublado.org/. 



absolute upper limit to the BEL flux. We notice that Rodriguez - 
Ardila, et al. (2000) found that on average, 50% of the flux of 
the total H/3 is due to emission from the NELR, and that the 
F([0 III]A5007)/ F(H/3) emitted in the NELR varies from 1 to 
5, which is much lower than our adopted value (~ 8.3). This 
result is sustained by the analysis of Contini, Rodriguez- Ardila, 
& Viegas (2003). However, Veron-Cetty, Veron, & Goncalves 
(2001) point out that the low F([0 III]A5007)/ F(H/3) values 
found by Rodriguez-Ardila, et al. (2000) for NLSyls are due 
to the fact that they modeled the broad Balmer component 
with a Gaussian rather than a Lorentzian. In the analysis of 
Veron-Cetty, Veron, & Goncalves (2001) the F([0 III]A5007)/ 
F(H/3) ratios span the range measured in BLSys, which 
makes our use of a value derived from a well-studied BLS1 
reasonable. Furthermore, Nagao, Murayama, & Taniguchi 
(2001) also show that F([0 I]A6300)/F([O III]A5007), 
F( [O III] A4363)/F( [O III] A5007) are indistinguishable in NLS 1 s 
and BLSls. Given this ambiguity, in the following analysis 
we will note where our assumptions for the deconvolution of 
O VIAA1032J038 and C IIIA977 affect the results. 

As a comparison, Column (4) reports the corresponding val- 
ues for a mean QSO spectrum, which we derived from Baldwin 
et al. (1995) by applying the reddening correction appropriate 
for Ark 564; the Ha/Lya ratio is derived from Osterbrock & 
Pogge 1985). Column (5) lists the BEL fluxes of the Syl.5 
NGC 5548, corrected for NEL contribution and Galactic red- 
dening (Korista & Goad 2000 and references therein). Column 
(6) lists the FWHM of the lines, as drawn from Papers V, II and 
III (C IIIA977 and O VI are the model BEL and NEL compo- 
nents, while the others are measured on the whole line profile). 
All errors are propagated in quadrature. Finally, Column (7) 
reports the references for Columns (2), (3), and (6). Table 3 
shows that Carbon in Ark 564 is at the lower end and Nitro- 
gen at the upper end of the mean QSO distribution, and some 
interesting differences can be found with respect to the Syl.5 
NGC 5548. Indeed, N V A 1240 is stronger in Ark 564 by a 
factor of - 2.3, while C IV A1550, C III] A1909, and Mg II 
A2800 are weaker in Ark 564 by a factor of - 4.4, 2.5 and 2.8, 
respectively. 

4.2. Input Continua 

We used the code Cloudy 18 (v94.00, Ferland 1996) to pre- 
dict the intensities of the lines produced by the BEL gas through 
photoionization modeling. Our choices of input continua for 
Cloudy are shown in Figure 4. In brief, 

1. The Cloudy "table agn" continuum, which is the 
Mathews & Ferland (1987) continuum modified with a 
sub-millimeter break at 10 /mi, so that the spectral in- 
dex is changed from -1 to -5/2 for frequencies below 
the millimeter break. 

2. The SEDA input continuum, which was created from 
points chosen from the spectral energy distribution pre- 
sented in §3 (the circles in Figure 4), as well as one 
extrapolated point (the empty square), i.e., where the 
optical/FUV and X-ray power law extrapolations meet. 
In particular, in the X-ray, we used continuum points 
from the power-law fit and added a black body com- 
ponent (the dot-dash line in Figure 3) of temperature 
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T = 1.8 x 10 b K and luminosity L bb = 2.48 x 10 J8 ergs 
s-'(§2.2). 

3. The more conservative SEDB input continuum, which 
only uses points from the observed spectral energy dis- 
tribution (circles only). 



lOOAim 10 M m l M m IQOOA O.lkeV IkeV lOkeV lOOkeV 




15 16 17 
log i/(Hz) 



20 



FIG. 4. — Input SEDs for Cloudy, normalized to the absorption-corrected 
rest-frame flux at 2 keV. The dotted line is the "table agn" model in Cloudy; 
the short-dashed line is the SEDA described in §3 (the circles are the observed 
points, the empty square the extrapolated point); the solid line is the conserva- 
tive SEDB described in §3. 
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FIG. 5. — Line intensity relative to Lyo of C III] A 1909 as a function of hy- 
drogen density (log n = 9, 1 1 , 12, 13) and input continuum ("table agn", SEDA, 
and SEDB discussed in §4.2). The horizontal dash-dot lines correspond to the 
observed values and their errors. 
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4.3. Physical Conditions of the Emission-Line Gas 

Figure 5 shows the line intensity of C III] A 1909 relative 
to Lya as a function of hydrogen density (log n (cm -3 ) = 
9,11,12,13) and input continuum ("table agn", SEDA, and 
SEDB discussed in §4.2). The horizontal lines correspond to 
the observed value of C III] A 1909 and its errors listed in Ta- 
ble 3. Solutions to U (as also shown in Figures 6 and 7) are 
double valued, but we prefer higher values based on line widths, 
as discussed below. As expected for a semiforbidden transi- 
tion, the C III] A 1909 line becomes collisionally suppressed as 
the density increases, arguing for an upper limit for the den- 
sity of log n £12. Analogously, Figure 6 shows the intensity of 
C IIIA977 and total N VA1240 which indicate logn > 9. Fur- 
thermore, the observed C III] A1909 implies logf/ ~ -3.3 or 
logf/ i=a -0.7 for \ogn = 9, both of which are much lower than 
the values required for O VI, for which we derive log U ~ -1 . 1 
or 0.8 19 (from a plot analogous to Figure 7b relative to logn = 
9), i.e., no consistent result can be found for low densities. 

Considering log n s=a 11 as a plausible estimate of the den- 
sity, we can investigate the value of the ionization parame- 
ter. This has been a difficult problem always, because with 
only optical/UV observations, multiple ionization states of a 
single element are not observed and so the values of U are 
highly model dependent. With the multiwavelength, multi- 
mission observations of Ark 564, we now have observations 
of both C III (with FUSE) and C IV (with HST), so we can 
actually measure the ionization parameter. As shown in Fig- 
ure 7a, the C IIIA977/C IVA1550 ratio constrains the ioniza- 
tion parameter to logf/ = [-2.88,-0.22] 20 . Studies of rever- 
beration mapping have shown that the BLR is stratified (Pe- 
terson & Wandel 1999), and so a single value of U cannot 
possibly correspond to the entire BLR, but the range of U 
determined above must be the dominant range of U where 
C IV emission is produced. The higher ionization lines, e.g. 
O VIAA1032, 1038 and He IIA1640 are likely to be produced 
closer in with higher ionization parameter. Indeed, as shown in 
Figure 7b, and Figure 6d, somewhat higher values of U are pre- 
ferred for O VI (log U = [-1 .52, -1 .36] or log U = [0.61 , 0.87] 21 ) 
and NV (logf/ = [-1.43,0.4]). 

Figure 7c shows the C IVA1550 intensity ratio and implies 
that either logf/ = [-2.94,-2.66] or logf/ = [0.18,0.34]. Fig- 
ure 7d shows the intensity ratio C IVA1240/N VA1550, for 
which we derive logf/ = [0.07,0.39]. The comparison with the 
BLS1 NGC 5548 shows that N V is stronger in Ark 564 by a 
factor of <~ 2.3, while C IV is weaker in Ark 564 by a factor of 
~ 4.4 (§4.1 and Table 3). Hence, the observed C IV/N V ratio 
may be roughly one order of magnitude smaller in Ark 564, and 
the observed limits in Figure 7d probably reflect an overabun- 
dance of N (or C depletion). 

To discriminate between the low-f/ (logf/ w — 1.5) and high- 
U (logf/ ss 0) solutions, we derive the distance of the BELR gas 

from Kbelr = (Q/4iTcnU) 1 ' 2 w 1.63 x lO^iQss/U) 1 ' 2 ^^ 2 
cm, where Q55 = (2/10 55 , assuming the above values of U and 
n, and the photon luminosity Q from §3. For logn =11 and 
logf/ w -1.5, the inferred distance of the BELR is = 
9.2 x 10 16 cm for g(SEDA), and ffffgf = 6.5 x 10 16 cm for 
g(SEDB). For logn = 11 and logf/ w 0, R s ™* = 1.63 x 10 16 



cm and R s B E E fl ■ 



1.15 x 10 16 cm. For a central mass of M - 



19 We find log U ~ -0.9 or 0.8 if we do not deconvolve BEL and NEL. 

20 logf/ = [-3.07,-0.36] if we do not deconvolve BEL and NEL. 

21 logf/ = [-1.44,-1.23] or logf/ = [0.64,0.79] if we do not deconvolve BEL and NEL. 
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FIG. 6.— Same as Figure 5 for C IIIA977 (top) and N VA1240 (bottom) for 
log n = 9 and 11. 
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FIG. 7.— Line intensity ratios of C IIIA977/ C IVA1550 (a), O VI/Lya 
(b), C IVA1550/Lya (c), and C IVA1550/N VA1240 (d), as function of input 
continuum (same line notations as in Figure 5). For all cases, log n = 11. 

8 x 10 6 M Q , the expected velocity dispersion (Peterson et al. 
2000) is V = (GM/R BLR f) - 5 , with /= 3/V2. This corresponds 



to 740-880 km s" 1 for the low-U solutions and 1750-2090 km 
s" 1 for the high-£7 solutions. Column (6) of Table 3 shows that 
FWHM(C IV)= 1934 km s _1 , FWHM(C IH)= 1920 km s~\ and 
FWHM(N V)= 2809 km s" 1 . Therefore, the comparison with 
the observed FWHMs favors the high-f/ solutions. Finally, Ta- 
ble 3 shows that the FWHM of O VI is larger than that of N V, 
which is larger than that of C IV, again consistent with the strat- 
ified BLR model. 

5. DISCUSSION 

A non-simultaneous optical, UV and X-ray spectral energy 
distribution of Ark 564 was presented by Comastri et al. (2001) 
who found that it peaks in the soft-X-ray band. Here we present 
a spectral energy distribution which is obtained from contempo- 
raneous data covering almost 5 decades in energy 22 . Simultane- 
ity is particularly important for NLSls, since, as a class, they 
are extremely variable in time, although Ark 564 has shown 
only weak variability in the optical/UV bands (Paper II-III). 

We report some relevant data from the spectral energy distri- 
bution in Table 2. These were derived using the simple parame- 
terization of the combination of an optical/FU V power law with 
spectral index a = 0.42 breaking at ~ 50 eV to ocasca = 1-538 
(SEDA in Figure 3), as well as the IRAS and IRTF data points 
(§2.1). A more conservative spectral energy distribution, in- 
stead, connects the FUV and the soft X-ray data with a simple 
power law (a = 1.08, the solid line in Figure 3), and is a com- 
bination of an optical/FUV power law with a = 0.42 breaking 
at 1000A to a = 1.081, then again breaking at w 0.8 keV to 
otASCA = 1-538 (SEDB in Figure 3). The ambiguity in the shape 
of the spectral energy distribution in the 900A-0.8keV region 
is rather unfortunate, since a considerable portion of the en- 
ergy of Ark 564 might be output in this range. Previous ROSAT 
(Brandt et al. 1994) and BeppoSAX (Comastri et al. 2001) data 
showed a flattening of the soft excess toward the lowest X-ray 
energies available, and the XMM spectrum obtained during the 
monitoring campaign of 2000 shows a definite curvature in the 
soft excess (Vignali et al. 2003). 

An interesting issue is how Ark 564 compares to other NLS 1 
galaxies and with BLS 1 in terms of its broad-band properties, 
as they can be quantified by spectral indices. Table 4 reports 
the intrinsic spectral indices calculated between different wave- 
length bands for Ark 564 and, as a comparison, the correspond- 
ing values for the NLS1 Ton SI 80 and BLS Is. These show 
that while the inter-band properties of Ton SI 80 are not signif- 
icantly different from the ones observed in BLS Is (Turner et 
al. 2002), this may not be the case for Ark 564. Table 4 indi- 
cates that the two NLS 1 s have steeper X-ray slopes than BLS 1 s, 
Ark 564 more so than Ton SI 80 (ct x = 1.57 and 1.44, respec- 
tively, compared to 0.91 for BLS Is), which is consistent with 
the general characteristics of NLSls. The optical/X-ray spec- 
tral index (a ox ), on the other hand, is lower in Ark 564 than in 
BLS Is and Ton S 180 with a ox = 1.11 for Ark 564 23 , 1.52 for 
Ton SI 80. This reflects the fact that Ark 564 is relatively more 
X-ray bright, or that the optical continuum is suppressed, com- 
pared to other AGNs. The other indices also reflect the X-ray 
brightness of Ark 564 when compared to BLS Is. We cannot 
exclude the possibility that intrinsic reddening in excess of the 



22 The Lick spectrum, which was obtained in 1980 is merely used to fill in a small gap in the optical data, and the 2001 FUSE spectrum shows a continuum slope 
consistent with the one obtained from the HST spectrum. This latter fact suggests that, although the flux level changed between 2000 and 2001, the overall shape of 
the optical/FUV spectral energy distribution did not change. 

23 The observed o ox is 0.94 in Ark 564. 
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FIG. 8. — Comparison of the SEDs of Ark 564 with the mean SED for radio- 
quiet quasars (RQ; Elvis et al. 1994) and LZ (Laor et al. 1997; Zheng et al. 
1997), one Seyfert 1 galaxy (NGC 4151, Kraemer et al. 2000), one Seyfert 1.5 
galaxy (NGC 5548, Kraemer et al. 1998, and the NLS1 Ton S180 (Turner et 
al. 2002). [A color version of this plot is available in the electronic edition of 
the Journal.] 

one considered here may be present that could be suppress- 
ing the optical, but we argue that it would require the contin- 
uum to be more reddened than the emission lines (e.g. He II 
A1640/A4686 ratio; Paper IV) show. If these differences, how- 
ever, are not due to intrinsic reddening, then they represent a 
true difference in the energy balance of this NLS 1 and that the 
X-ray emission is therefore consistent with the predictions of 
the slim disk models, indicating high accretion rates. 

We estimated the luminosity in different energy bands us- 
ing our parameterizations of the spectral energy distribution de- 
scribed above (SEDA and SEDB). Though with a ~ 20% sys- 
tematic error in the 0.3-0.6 keV band flux (C. Vignali 2003, 
private communication), the XMM spectrum provides a way to 
extend the energy range below the ~ 0.8 keV limit of the ASCA 
spectrum. Therefore, we also estimated the luminosity for the 
best-fit model to this spectrum. Table 5 reports these luminosi- 
ties, and indicates that more energy is emitted in the 10-100 eV 
band than in any other decade, constituting roughly half of the 
emitted energy in the opt/X-ray ranges. This implies that the 
primary spectral component peaks in the extreme UV-soft X- 
ray band and this is generally consistent with disk-corona mod- 
els (Haardt & Maraschi 1991). The integrated luminosity be- 
tween 10~ 5 and lOkeV gives a lower limit on the bolometric 
luminosity of L\, \ ^ 10 45 ergs s" 1 . 

Growing evidence has been gathered that NLS Is are super- 
Eddington accretors. Collin & Hure (2001) and Collin et al. 
(2002) studied accretion rates in a sample of AGNs using the 
Kaspi et al. (2000) relationship, and found that not only half 
of their sample accretes close to or above the Eddington rate, 



but also that the largest Eddington ratios are found in NLS Is. 
Indeed, they showed that NLS 1 s are at the extreme of a well- 
defined sequence relating the Eddington ratio to the line widths. 
Furthermore, Wang & Netzer (2003) show that a model of ex- 
treme slim disk (which is responsible for the soft X-ray excess, 
or hump, seen in most NLS Is) and a hot corona (contributing 
to the hard X-ray emission) can also naturally explain the X-ray 
spectral variability characteristics observed in Ark 564 such as 
simultaneous variations of the soft hump and the hard X-ray 
without a significant time lag (Turner et al. 2001). Through a 
comparison of the X-ray variability of Ark 564 and the BLS 1 
NGC 3516, Pounds et al. (2001) estimate that the mass of the 
central black hole in Ark 564 is ~ 10 7 M©, implying an accre- 
tion rate in the range m w 0.2-1. With a bolometric luminos- 
ity in the order of 10 45 ergs s" 1 and an Eddington luminosity 
^Edd ~ 10 45 ergs s" 1 , we also infer that m ss 1. Wang & Netzer 
(2003) also provide a means of estimating the black hole mass 
(if the accretion is super-critical) which is independent on the 
accretion rate itself. For Ark 564, we obtain M w 2 x 10 6 M Q , 
which is within a factor of a few from the Collier et al. (2001) 
estimate. 

Figure 8 compares the SED of Ark 564 with the mean SED 
for radio-quiet quasars (Elvis et al. 1994), the LZ SED (Laor et 
al. 1997; Zheng et al. 1997), the Seyfert 1.5 galaxies NGC 5548 
(Kraemer et al. 1 998) and NGC 4151 (Kraemer et al. 2000), and 
the NLS1 Ton SI 80 (Turner et al. 2002). There are significant 
differences in the intrinsic shape of the SED across the AGN 
population (see, also, Turner et al. 2002), the most evident be- 
ing the energy of the peak and the presence (or lack of) of the 
big blue bump (BBB), the signature of the emission from the 
accretion disk. The radial dependence of the temperature for an 
optically thick, geometrically thin accretion disk (Shakura & 
Sunyaev 1973) is, T(R) ~ 6.3 x 10 5 (m) 1 / 4 M~ 1/4 (R/R S )~ 3 / 4 K 
(Peterson et al. 2000), where M% is the mass in units of 10 8 M Q , 
R is the radius, and Rs is the Schwarzschild radius. For Ark 564, 
using m w 1, as we derived above, and M«8x 10 6 M Q (Pa- 
per II), the peak temperature is ~ 125 eV; this is within a factor 
of 3 from the peak of our less conservative parameterization, 
SEDA (defined extrapolating the optical power-law continuum 
to meet the extrapolation of the X-ray power law) which peaks 
at 50 eV. The true SED probably peaks somewhere between 
these two values. For comparison, the NLS Is RE J 1034+396 
(Puchnarewicz et al. 2001) and Ton SI 80 (Turner et al. 2002) 
peak at rj 250 eV and £ 100 eV, respectively. Therefore, even 
among the NLS Is, differences in the shape of the SED are ob- 
served. However, in none of these NLS Is there is an indication 
of the presence of optical/UV BBB, and a strong soft X-ray ex- 
cess is seen, instead. In this light, Ark 564, is also consistent 
with the paradigm that the accretion disk is so hot in NLS Is 
that the BBB is shifted in the EUV-soft X-rays. 

We also note that Ark 564 is rather FIR bright (with respect to 
the optical), compared to the sample of radio-quiet quasars and 
the LZ sample, as can be seen in Figure 9, where the SEDs have 
been normalized to match their optical/UV slope (as opposed to 
the 2 keV flux in Figure 8). Indeed, if we use the definition of 
the IR flux as a function of the IRAS fluxes given by Sanders 
& Mirabel (1996) 24 , we obtain that L(8- 1000/t m) « 1O H L , 
which makes Ark 564 a luminous IR galaxy. The shape of the 
IR/optical SED of Ark 564 also resembles the shapes observed 
in the IRAS Bright Galaxy Survey (see Figure 2 in Sanders & 
Mirabel 1996, for intermediate values of f^y). Using SEDA, we 



F(8-1000^m) = 1.8 X 10~ 14 {13.48/i 2 + 5.16/25 + 2.58/ 6 o+/ioo} W mT 2 , where fn,f2S,fm, and/ioo are the IRAS flux densities in Jy at 12, 25, 60 and 100 /im. 
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FIG. 9. — Comparison of the Ark 564 SED with the mean SED for radio- 
quiet quasars (RQ; Elvis et al. 1994) , LZ (Laor et al. 1997; Zheng et al. 1997), 
and the "table agn" model in Cloudy (AGN). The SEDs are normalized so that 
they match in the optical/UV part. 

estimate that the FIR luminosity is L FIR rs 1.3 x 10 44 ergs s" 1 , 
i.e., ~ 10% of the total luminosity and <~ 20% of the combined 
optical/UV/X-ray luminosity. Crenshaw et al. (2002) noted that 
the associated warm UV absorber is lukewarm and dusty. The 
IR emission observed in this object could then be thermal emis- 
sion from the dust grains embedded in the absorber as they are 
heated by the strong UV/EUV continuum. Given the IR bright- 
ness in this object and in many NLSls (Moran, Halpern, & 
Helfand 1996), it is not unlikely that a contribution might be 
coming from the host galaxy, in the form of a nuclear starburst 
(Mathur 2000). Crenshaw et al. (2002), however, did not de- 
tect any extended emission in the two-dimensional STIS spec- 
tral images that could be due to a nuclear starburst. We com- 
pared the FUSE spectrum, which was taken through a much 
larger aperture, and hence is more likely to show stellar ab- 
sorption features, with a FUSE 'template' spectrum of the star- 
burst galaxy NGC 7496; using the constraints from the Ly7- 
C III] A977 profiles, we obtained a rough upper limit of the star- 
burst contribution at 1000A of 50%. A further constraint on 
the IR starburst contribution would probably come from detec- 
tion and measurements of the 3.3-3.4/im PAH emission fea- 
tures, which have been found in starburst galaxies, luminous 
IR galaxies, and obscured AGNs (Moorwood 1986; Imanishi 
2002), and which have been successfully detected in the NLS 1 
NGC 4051 (Rodriguez- Ardila & Viegas 2003). 

Using our SEDA and SEDB as inputs to Cloudy we pre- 
dicted the intensity of the strongest lines in the FUV/UV spec- 
trum and compared them with the observed values in order 
to constrain the physical parameters of the line-emitting gas, 
namely, the density and ionization parameter. Figure 6 and 7 
show that SEDA and SEDB, because of their strong EUV to 
soft X-ray flux, for a given observed line ratio predict values 
of ionization parameter which are lower than those with stan- 



dard AGN continuum. From C III] A 1909, C III A977, and N V 
A 1240, we infer that logn w 11. Two classes of solutions for 
U are consistent with this density value, one with low U val- 
ues (logf/ s=a -1.5) and one with high U values (logf/ « 0). 
We discarded the low-f/ class (§4.3) on the basis that the pre- 
dicted widths of the lines, derived from the velocity dispersion 
V = (GM/R BLR f)°\ of 740-880 km s" 1 are too small with re- 
spect to the observed ones (« 2000 km s" 1 ; Paper II). As ex- 
pected, we find that the BLR is stratified around log (/kO, with 
higher ionization lines originating from regions with higher U . 

Column (6) of Table 3 shows that FWHM(Lya)= 2114 km 
s" 1 , FWHM(N V)= 2809 km s" 1 , and FWHM(C IV)= 1934 km 
s" 1 , which indicate that the radii of the La, N V, and C IV 
broad-line emitting regions are R\£r ~ 4.3 lt-days, R^r ~ 2.5 
lt-days, and R™ R ~ 5.2 lt-days. Using the findings of pre- 
vious monitoring programs on Seyfert Is (Netzer & Peterson 
1997), we can estimate the size of the H/3-emitting region 
from R^ R w 0.5R^ R , and R^ R w 0.2R^ R . For NGC 5548, 
furthermore, R™ R w 0-5R^[ R (Peterson 1993). We can con- 
clude that /?blr ~ 10 ±2 lt-days, which is consistent with the 
/?gL R -luminosity relationships of Kaspi et al. (2000) and Pe- 
terson et al. (2000), when we assume a luminosity ALa(5100 
A)w 3.2 x 10 43 ergs s" 1 (Table 2). This indicates that the BLR 
radius of this NLS 1 is consistent with the distribution of BLR 
radius in BLSls, and that the narrowness of the emission lines 
is not due to the BLR being relatively further away from the 
central mass than in BLSls of comparable luminosity. 

Table 3 shows that some interesting differences in line ratios 
can be found with respect to the Syl.5 NGC 5548. Indeed, N V 
A1240 is stronger in Ark 564 by a factor of <~ 2.3, while C IV 
A 1550, C III] A 1909, and Mg II A2800 are weaker in Ark 564 by 
a factor of ~ 4.4, 2.5 and 2.8, respectively. While all line ratios 
in Ark 564 are statistically consistent with the ones measured 
for a mean QSO (given the large uncertainties on our measure- 
ments), Carbon lines are at the lower end and Nitrogen at the 
upper end of the QSO distribution, confirming this trend for 
weak Carbon and strong Nitrogen in Ark 564. Furthermore, 
C III] A977 would indicate that -3.13 < logf/ < -0.05 (Fig- 
ure 6), and for this range, the observed N V/C IV ratio is larger 
than the model predictions by a factor of ~ 8. This may imply 
super-solar metallicity in this NLS1 as suggested by Mathur 
(2000) and is consistent with the finding of Shemmer & Netzer 
(2002) that NLSls have higher metallicities than BL AGNs for 
a given luminosity. 

An interesting question is how sensitive the emission lines 
are to the true shape of the ionizing continuum. Given the dif- 
ficulty related to deblending the BEL and NEL components of 
the lines and consequent large errors on the observed line ra- 
tios, the present emission line data do not allow us to discrimi- 
nate between SEDA and SEDB. These SEDs differ in the range 
1000A-750eV (corresponding to the gap in the data between 
the FUSE and ASCA spectra), with the maximum difference at 
around 50 eV. This difference should show the most for the high 
ionization lines of C IV, N V and O VI. However, the predicted 
strength of their emission lines is very similar for the two SEDs, 
for a wide range of ionization parameters of interest (Figures 6 
and 7). Perhaps, this is the reason why the emission line spectra 
of most AGNs look so very similar, over a wide range of lumi- 
nosities. This demonstrates how unsuitable emission lines are 
as diagnostics for the underlying SEDs. 

The absorption lines, on the other hand, are sensitive to the 
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U 



FIG. 10. — Top panel shows the ratio of the column densities relative to 
SEDA and SEDB as a function of U for C IV, N V, O VI, and O VII. Bottom 
panel shows the ratio of the column densities of O VII to C IV as a function 
of U for solar abundances, for SEDA and SEDB. 

input SED (Mathur et al. 1994). Column densities of differ- 
ent ions can be inferred from the the fractional abundances 
/i on , the total column density Nh, and the assumed abundances 
Nx, Ni 0!1 = NnNxNion, through photoionization calculations (as- 
suming an SED). Figure 10 (Top) shows the fractional col- 
umn densities for the same ion calculated with the different 
SEDs [ATciv(SEDA) /AT C iv(SEDB), AW(SEDA) /AW(SEDB), 
AWi(SEDA) /W vi(SEDB), and Afovn(SEDA) /AWSEDB)], 
and illustrates how the predicted column densities of high ion- 
ization lines differ, with O VI showing the largest variations, 
with over two orders of magnitude difference. Figure 10 (Bot- 
tom) compares the column density for O VII and C IV calcu- 
lated with the same SED and solar abundances [TVovii(SEDA) 
/ JVciv(SEDA) and AWn(SEDB) / Afciv(SEDB)], which turn 
out to be good diagnostics even at lower ionization parameters. 
Thus, in principle, absorption line studies offer a far more pow- 
erful tool to determine the ionizing continuum of AGNs, com- 
pared to the emission lines. This also underscores the power of 
multiwavelength observations, as the maximum discriminator 
comes from comparing the lower- and higher-ionization lines. 



6. SUMMARY 

We presented the intrinsic spectral energy distribution of 
Ark 564, constructed with quasi-simultaneous data obtained 
during 2000 and 2001. We compared this SED with that of 
Ton SI 80 and with those obtained for Broad-Line Seyfert Is 
to infer how the relative accretion rates vary among the Seyfert 
1 population. The peak of the SED is not well constrained; 
however, in our parameterization most of the energy of this ob- 
ject is emitted in the 10-100eV regime, and constitutes roughly 
half of the emitted energy in the optical/X-ray ranges. This is 
consistent with a primary spectral component peaking in the 
soft X-ray band, therefore with the predictions of the slim disk 
models, hence high accretion rates. Indeed, we estimate that 
m s=s 1. 

We constrained the mean physical conditions in the BELR of 
this AGN, by examining the emission lines observed in its spec- 
trum, and deduced the physical properties of the line-emitting 
gas through photoionization modeling. We concluded that the 
line-emitting gas is characterized by log n w 11 and log U ~0, 
and is stratified around log U ~0. Our estimate of the radius 
of the H/3 emitting region R^lr ~ 10 ± 2 lt-days, is consistent 
with the /? B L R -luminosity relationships of Kaspi et al. (2000) 
and Peterson et al. (2000). This indicates that the narrowness of 
the emission lines is not due to the BLR being relatively further 
away from the central mass than in BLSls of comparable lumi- 
nosity. We also find evidence for super-solar metallicity in this 
NLS1, based on the low C IV/N V observed line ratio. While 
the emission lines turn out to be unsuitable as diagnostics for 
the underlying SEDs, we showed that absorption line studies 
offer a far more powerful tool to determine the ionizing contin- 
uum of AGNs, especially if comparing the lower- and higher- 
ionization lines. This underlines the power of multiwavelength 
observations. 
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Table 1 

Observing Log for Arakelian 564 



Observatory 


Instrument 


UT Dates 


Wavelength/Energy 3 


Notes 


References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


ASCA 




2000 Jun 1-M6 


0.75-9.76 keV 


2.98 Ms, continuous 13 


1 


XMM-Newton 




2000 Jun 17 


0.3-8 keV 




2 


FUSE 




2001 Jun 29-30 


1000-1 175 A 


63ks;30"x 30"(LWRS) 


3 


HST 


STIS/G140L 


2000 May 9-Jul 8 


1175-1711 A 


554304 s; 52"x0."5 


4,5 


HST 


STIS/G230L 


2000 May 9-Jul 8 


1711-3143A 


24216 s; 52"x0."5 


4,5 


Lick 




1980 


3170-4160A 




6 


Wise 


FOSC 


1998 Nov-2001 Jan 


4 160-7790 A 




7 


IRTF 


SPEX 


2000 Oct 11, 13 


8200-24000 A c 


-30min, 15"x0"8 


8 


IRAS 






12, 25,60, 100 yum 




9 



a Observed-frame wavelength/energy bands utilized. 

b Except for gaps due to Earth occultation and passage of the spacecraft through the SAA. 
c Only few continuum points were used for this work. 

References. — (1) Turner et al. 2001. (2) Vignali et al. 2003. (3) Romano et al. 2002. (4) Collier et al. 2001. (5) 
Crenshaw et al. 2002. (6) D. E. Osterbrock 2002, private communication; (7) Shemmer et al. 2001. (8) Rodriguez- Ardila 
et al. 2002b and references therein. (9) Moshir et al. 1990. 
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Table 2 



Data from the Spectral Energy Distribution 



Rest Wavelength 


ul ' j f Observed^ 


vL v (Intrinsic, SEDA) 


/Energy 


(x 10 43 ergs s -1 ) 


(x 10 43 ergs s -1 ) 


(1) 


(2) 


(3) 


97.59 Aim a 


3.98 


4.57 


58.56 /zm a 


4.82 


5.50 


24.40 /xm a 


7.91 


8.91 


11.71 /zm a 


8.73 


9.77 


2.4 ^m b 


3.3 


3.5 


i .0 /im 


Z.z 


1 A 


1.14^m b 


2.1 


2.4 


1 /^m c 


2.045 


2.165 


9850,A b 


2.0 


2.3 


7000 A 


1.851 


2.660 


5500 A 


1.73 


3.057 


5 100 A 


1.694 


3.193 


3000 A 


1.461 


4.338 


2500 A 


1.388 


4.819 


1000 A 


1.075 


8.178 


0.046 keV d 




17.5 


0.25 keV d 




7.195 


0.78 keV 


2.895 


5.240 


IkeV 


3.046 


4.668 


2keV 


2.060 


2.419 


lOkeV 


0.942 


0.970 



Note. — The intrinsic optical/UV/X-ray data are from the 
reddening-corrected, rest-frame SEDA (§3). SEDB is SEDA 
with the exclusion of the point at ~ 50 eV. We adopt Hq = 75 
km s" 1 Mpc" 1 , qo = 0.5. 

*IRAS data points (§2.1). 
b IRTF data points (§2.1). 

c Extrapolated value from the optical/FUV power law (spectral 
index a = 0.42 ±0.01, §3). 

d Extrapolated value from the ASCA power law (a AS cA = 1.538 ± 
0.005, §2.2), peak of SEDA (§3). 
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Table 3 

Arakelian 564 Emission-Line Characteristics 



Line 


F/F(Lya) 


F/F(Lya) 


F/F(Lya) 


F/F(Lya) 


FWHM 


References 




DCI _i_MCT a 


TJT7T a 




(Km s ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


C III A977 


0.053±0.017 


0.031±0.015 






4000,1100 


1,2 


OVI A1032 


0.284±0.032 


0.200±0.052 




0.036 


4000,1100 


1,2 


OVI A1038 


0.141±0.048 


0.100±0.026 




0.018 


4000,1100 


1,2 


VI A1032,1038 


0.425±0.141 


0.300±0.076 


0.1-0.3 e 


0.054 


4000,1100 


1,2 


La A1216 


LOO 


1.00 


1.00 


1.00 


2114 


3,4 


N V A1240 


0.275±0.118 




0.09-0.26 


0.119 


2809 


3,4 


CIV A1550 


0.215±0.087 




0.28-0.42 


0.937 


1934 


3,4 


Hell A 1640 


0.101 ±0.040 




0.06-0. 13 f 


0.143 f 


1195,1831 


3,4 


III] A 1663 


0.038±0.015 










3,4 


CIII] A 1909 


0.068±0.026 g 




0.09-0. 19 h 


0.171S 


1920 


3,4 


Mg II A2800 


0.066±0.008 




0.06-0.13 


0.188 


1659 


3,5 


H/3 A4861 


0.063±0.019 




0.02-0.05 




700 


3 


Ha A6563 


0.240±0.080 




0.05-0.09' 






3 



a Reddening-corrected flux relative to Lya derived from Crenshaw et al. (2002). The lines are corrected 
using E(B-V) = 0.14 ± 0.04 mag and Ark 564 reddening curve from Crenshaw et al. (2002) plus E(B-V) = 
0.03 mag and Galactic curve (§2.3). 

b Data from the mean observed QSO spectrum in Baldwin et al. (1995) and references therein, corrected for 
the reddening appropriate for Ark 564. 

C BEL fluxes, corrected for NEL contribution and Galactic reddening of the Syl.5 NGC 5548; derived from 
Korista & Goad (2000) and references therein. The O VI line ratios are derived from the BEL values in 
Brotherton et al. (2002), then corrected for Galactic reddening (E(B - V) = 0.03 mag, and extinction law of 
Cardelli, Clayton, & Mathis 1989). 

d Model FWHM of C IIIA977 and O VI relative to BEL and NEL components (Paper V), separately; the 
others are measured on the whole line profile (Paper II). The He II values are relative to the G140L and G230L 
mean spectrum, respectively (Paper II). 

Total O VI+Ly/3 flux. This compares with 0.434 ±0.141 for Ark 564. 

f Total He II+O III]A1666 flux. This compares with 0.139 ± 0.043 for Ark 564. 

Total C III] A1909+Si III] A 1892 flux. 

h Total C III] A1909+Si III] A1892+A1 IIIA1990 flux. 

'Based on the range of values of Ha/H/3 (3.97-6.64) from the NLS1 sample of Osterbrock & Pogge (1985). 



References. — 



(1) This work. (2) Paper V. (3) Paper IV. (4) Paper II. (5) Paper III. 
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Table 4 
Spectral Indices 



Index 


Definition 


Ark J04 


Ion MoU 




References 




(2) 


(3) 


(4) 


(5) 


(6) 


1 00 fim— 1 2 fim 


-1.086 log(Fi2 M m/^100Mm) 


0.64 








^*12 ^im— 2.4 /^m 


-1.431 l0g(F2Afim/Fi2fj, m ) 


1.64 








tt2.4pm-1.6/jm 


-5.679 \og(F L6t , m /F 2 A flm ) 


1.94 








ttl.e^im-l 


-4.900 logtF^/F.e^) 


1.18 b 








Ol2^m— 1 ^im 


-0.927 log(F 1(Um /F 12 ^ m ) 


1.60 b 








a 3000-1000 (<* U v) 


-2.096 log(F 100 o/F 3 ooo) 


0.42 


0.66 


1.25 


1 


aFUSE-ASCA C 




1.08 






1 


0-5500-0.25 keV 


-0.489 log(F .25kev/F 5 50o) 


0.82 b 


1.12 


0.73 


2 


0:5500-1 keV (^ox-hard) 


-0.378 log(F lkeV /F 550 o) 


0.93 


1.38 


1.13 


3 


<2l/im-2keV («ix) 


-0.312 log(F 2keV /F lM ) 


1.01 b 


1.35 


1.14-2.16 


4 


O2500-2keV (c*ox) 


-0.384 log(F 2ke v/F 25 oo) 


1.11 


1.52 


1.46^$, 1.21±0.02 


5,6 


o x 


-1.431 log(Fiokev/F 2k ev) 


1.57 


1.44 


0.91 


7 



Note. — For spectral indices relative to the Ark 564 IRAS points we used the reddening-corrected, redshift- 
corrected fluxes and wavelengths. 

intrinsic, i.e, reddening-corrected, redshift-corrected (SEDA). 
b Based on extrapolated value (§3, Figure 3). 

c Simple power law connecting the high energy end of the FUSE spectrum and the low energy end of the ASCA 
spectrum (§3, Figure 3). 

References. — (1) Cheng, Gaskell & Koratkar 1991; index in the 2200-1200 A band, based on the BLS1 sub- 
sample. (2) Turner et al. 1999a. (3) Grupe et al. 1998. (4) Lawrence et al. 1997. (5) Zamorani et al. 1981. (6) 
Puchnarewicz et al. 1996. (7) Nandra et al. 1997b. 



Table 5 
Luminosities 

Energy Range L (Observed) L (SEDA) a L (SEDB) a L (XMM) a 

(keV) (xlO^ergss" 1 ) (x lO^ergs s" 1 ) (x 10 44 ergs s" 1 ) (x 10 44 ergs s" 1 ) 



(1) 


(2) 


(3) 


(4) 


(5) 


10- 5 -10" 4 




2.1192 


2.1192 


2.1192 


10- 4 -10" 3 




1.8108 


1.8108 


1.8108 


10 _3 -0.01 




1.2464 


1.2464 


1.2464 


0.01-0.1 




3.5175 


2.2589 


4.1430 


0.1-1 


0.2201 


1.6314 


1.8080 


2.7716 


1-10 


0.4086 


0.4485 


0.4485 


0.5285 


10" 5 -10 




10.7738 


9.6918 


12.6202 



a Reddening-corrected, rest-frame luminosities. 

Note. — The SEDA luminosities have been calculated using power-law parameterization of 
the SED with the spectral indices reported in Table 4. The SEDB luminosities refer to the more 
conservative parameterization described in §3 (Figure 3). The XMM luminosities make use of 
the XMM spectrum in the 0.05-10keV band. 



